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Recently, there have been considerable interests in theoretical studies of Kondo insulators 
and heavy fermion systems, the class of rare earth compounds [p]-p^. The mechanism for 



these systems to become insulating is due to the interaction between the conducting electrons 
and the array of localized impurity spins. On a one dimensional lattice, when the electrons 
hop between nearest-neighboring sites, the metal-insulator phase transition in the system at 
half-filling described by the Kondo lattice model has been investigated with various methods 



||q JTT|JTq j7| JT7|JT6|] . It was found that the simple Fermi liquid fixed point of the system (J = 0) 
is unstable against infinitesimally small interaction. Arbitrarily small interaction would drive 
the system from conducting phase to an insulating phase. 

More recently, Carruzzo and Yu develop a novel self-consistent variational approach, to 
investigate again whether an external magnetic field could induce a metal-insulator phase 
transition in the Kondo insulators 0. It is found that the insulating gap never vanishes 
for any value of magnetic field, suggesting no field-induced metal-insulator transition, in 



contrast to the predictions by the slave boson mean field theory [p!^-p!^,p[] . 

In the self-consistent variational theory, two scattering processes between the electrons 
and the impurity spins are kept in consideration: the process of no-spin-fiip and zero mo- 
mentum transfer (Si), and that of spin-flip and vr momentum transfer (5*2). In absence 
of external magnetic field, the approach has indicated that the insulating phase is stable 
even for arbitrarily small interaction parameter, that is, Jc = O'^, with a small but finite 
energy gap decaying exponentially, consistent with the previous results obtained through 
other methods 0. In the following, we use their approach to study the stability of the 
insulating phase of the ID chiral Kondo lattice model introduced recently by us. We find 
that at half-filling, the insulating phase is stable for any nonzero coupling constant. The 
insulating gap is computed explicitly for this system. 

The ID chiral Kondo lattice model consists of conducting electrons moving in one direc- 
tion on a close chain. At each lattice site, the conducting electrons interact with a localized 
impurity spin through spin exchange interaction. The system is described in terms of the 
following Hamiltonian: 



= E E <k)clc,„ + JE4^Q, ■ - /.E[4q„^ + 4(0], (1) 

where the conduction band spectrum is e{k) = —tk, with — 7r(l — < k < 7i{l—/L), in the 
momentum space. J is the couphng constant between the local impurity moments and the 
conducting electrons, and h is the external magnetic field. The local moments are described 
by the spin 1/2 operators, that is, [SJ{i), S^{i)] = iS^{i) (plus two other commutation 
relations obtained by the cyclic permutations of x,y,z), with the relation Sj{i) =3/4, for 
all the sites i = 1,2, ■ ■ ■ , L. In the following, for simplicity, we always restrict ourselves to 
the parameter range t > and antiferromagnetic interaction J > 0. The dispersion relation 
of the conducting band e{k) = —\t\k yields the electron group velocity v = de{k)/dk = —\t\. 
Below, we only consider the thermodynamic limit L — oo. 

In the extreme limit J = +oo, away from half-filling, we have found that the wavefunc- 
tions of the system are the RVB-type Jastrow product wavefunctions The interesting 
aspect is that in this limit various correlations of the system can be computed exactly in 
compact form. Moreover in this limit, the system obviously is a conductor away from half- 
filling, while it is an insulator at half-filling. For very large but finite J and at half-filling, 
each impurity spin will also attempt to form a singlet with one conduction electron at each 
site, to lower the energy of the system as much as possible. To transfer one electron with 
the long range hopping matrix element (corresponding to the chiral band) from one site to 
another would break two singlets, giving rise to a charge gap of the order 0{J), making 
the system an insulator. One interesting issue is the stability of the insulating phase. In 
the following, we'll investigate the stability of the insulating phase of the system when one 
changes J, and we'll restrict ourselves to the case of half-filling. 

Following Popov and Fedotov [^], one first introduces fermion operators {/|q,; fja},j G 
{1, 2, ■ ■ ■ , L}, a G {+, — } = {t, i} to describe the impurity spins, together with an additional 
term —-^J^jafjafja to take into account the constrain that every site is occupied by one 
impurity. The partition function of the system thus takes the following form: 

Z = J Dc^DcDpDf e-/o''^-etiE.=±4M9^--W+EtiE„.±/LM9./«.W+^^W]^ (2) 
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where c and / are Grassman variables, and H{t) is given in momentum space by 

H{r)= E E Hk) - ah/2]ci{k, T)c^{k, r) + E E " «V2]/I(fc, r)Uk, r) + 

k a=± k a=± '^P 

+ ^ E 5(^1 - ^2 + ^3 - /^4)/l(/i;i,T)//3(fc2,1")4(/i;3,T")c5(fc4,1")o'a/3 ■ (J^5. 

fci,---,fc4;«,/3,7,(5 

(3) 

Furthermore, introducing "right" and "left" operators for electrons and impurities, i.e. 

, , , = > 0, a = 

0a,a(A;,r) = <( (4) 

4>L;a{k, t) = 4)a{k — TT, r) A; > 0, a = 1 

with a G {R,L} = {0, 1}, and (f)a,a standing either for Ca,a or for fa^a, one can write H{t) 
as follows: 

^(^) = E E E [<k - an) - ah/2]cijk, T)ca,c.{k, t) + 
+ E E E[-^-«V2]/i,JA;,r)A,«(A:,r) + 
+ E E E fU^i^^)C%hAk2,r\ (5) 

A:i,fc2>0 a,6=0,l «,/3=± 

where 

^S = 77 E E E 4,7(^3, r)Cd,5(fc4,r) [<J«/3 ■ X 

X(5(A;i — k2 + k^ — ki — {a — h + c — d)7r), (6) 

with c?a/3 ■ c?^5 = a'-y6a,f3S^^s + 2(5q,,_/35^,_5, and we have neglected the k = terms. 
Following Ref. 0, we shall keep only the scattering processes: 

Si : ki = k2, k^, = k4,a = b,c = d,a = = S, 

82- ki = k2, ks = k^, a = b + 1, c = d + 1 {mod2), a = — /3, 7 = —6. (7) 

The two scattering processes correspond to the following effective interaction: 5*1 : 
^Ei,aa/iIa/i,J[Ej,7 7c],7Ci.7]''^2 : ^ [Ei,a (-l)7ila/i,-a] [Ei,7 (" l)^'4^Cj- . Therefore, C 
is the 4x4 matrix given in this approximation by 



c 



^ab Qab 

with matrix elements given by 

J 



C — \.(Jx — C_(_+l^ 



(8) 



fc>0 6=0,1 7=± 

= WrJ2J2 c-b,-ik, r)cb^+{k, r), 
^-^ k>Ob=± 

= WtYY r)cb,-{k, t). 

k>Ob=± 

Introducing 4- vector /''' = (/q^, flp /q^, /||), the interaction term is thus given by 

k>0 

We can then rewrite the partition function in terms of following functional integral 



(9) 



(10) 



Z^J Dc^DcDfDf exp{-[5„e[c^ c] + S^AfK f] + Sf,[c\ c, /t, /]]}, (11) 



where the various actions are given by 

Soc = ^ Efc>o Eu, c\k, uj)Doc{k, uj) 
[DqY^U = 5a,b5a,i3['^^ - "^/2 + e{k - an)] 



Sof = 4Efe>oE<^ f{k,^)Fof{k, 



UJ 



(12) 



iFo)^p = Sa,bSaA^uj-^-ah/2] 



and (j){k, t) = _ (2n+i)n (p{k, u;)e*'^'^. Keeping only the terms with Ui = UJ2 { thus = 004^ 



as considered in Ref. [01, we have 

Sfc = ^2 Efc>oEc. fKk,uj)Cf{k,uj) 



C 



C'h — <^x 



(13) 



with matrix elements: 



J 



^■^(^ k>0 6,7 
J 



or nY.Y.Y.^-b,~ik,uj)cb,+ ik, 

^^P k>0 ^ b 
^-^P fc>0 b 



00 



(14) 



For this chiral Kondo lattice model, we follow the self-consistent variational approach of 
Carruzzo and Yu. To find the quasi-particle spectrum of the electron degrees of freedom, 
one first integrates out the impurity degrees of freedom, to obtain the effective action of the 
electrons. When integrating out the Popov-Fedotov fermions for the impurity spins, one 
replaces the bare propagator with its dressed version, taking into account the effects of the 
conduction electrons on the impurity spins themselves. After some algebraic manipulation, 
we obtain the effective functional integral for the electron degrees of freedom: 



Z = y Dc^Dc exp —Seff[c\c], 



(15) 



with the effective action 



Seff = 7^ E E ^Wk, uj)c{k, u), (16) 

where using the same approximation of Ref. ||^, the corresponding matrix D is found to be 
(iu-h + e{k) \ 



D 








\ 4A 



iu — h + e{k — tt) 

iuj + h + e{k) 
iuj + h + e{k ~ n) J 








(17) 



where e{k) = —kt, Co and are the mean field parameters 



Co = lim 

/3-+00 



< C++ > _ < c+_ > 



Ctt = lim 



and 



J ,h 



X = [{h/2 - Cof + ClY'\ h = h/2- -(- - Co) 



(18) 



(19) 



Finding the inverse of the matrix D, with the analytical continuation iu —>■ w, the poles of 
the propagators yield the quasiparticle dispersion relations. We have found that there also 
exist four bands Ei{k) , E2{k) , E^lk) and E^lk), given below: 

E^{k) = kt- nt/2 - [{nt/2 + \h\f + {JCjAXfY^^ 
E2{k) = kt- 7rt/2 - [{7rt/2 - \h\y + (JC^/4A)Y/^ 
E^ik) = kt- nt/2 + [int/2 - \h\Y + (JCjAXf]'/^, 

E^{k) = kt- 7rt/2 + [(7rt/2 + \h\f + {JCjAXfY'\ (20) 



Let us only consider the case of no external magnetic field h = 0. For sufficiently large 
coupling constant J, when the system is in the insulating phase, the upper two bands E^{k) 
and Ei{k) are positive, the lower two bands Ei{k) and E2{k) are negative. The energy gap 
of the system is given by 

A = E,{0)-E2i7T). (21) 

The self-consistent equations for the mean field parameters Cq and C.,^ are given by: 

J'^ dk I 1 
~ 16 io 2^^^ ^ 

4 JO ivr £/4 £/3 

where A = {Cl + C2)i/2^ = JCo/(4A). If the quasi-particle energy gap A is greater than 
zero, the insulating phase is stable. However, as the coupling constant J is small enough, 
the energy gap approaches to zero, and the system becomes conducting at this critical point. 
In order for the theory to be self-consistent, the definitions of the mean field parameters 



C,r and Co must satisfy the condition \C-„\ < J/4, |Co| < J/4. The self-consistent Eq. ( p2D 
can be solved for Cq and C,r- One can find that Co = 0. The value of is determined by 
the following equation: 

dk 1 



iai= (J7i6) 



71 kt- nt/2 + [(7rt/2)2 + (J/2)2]i/2 



^/(l^l). (23) 



where the function f{x) = a;ln l±\l±^lJ—_ Since < /(x) < 1 for any a; > 0, one sees that 
ICttI < I, a physically acceptable solution. For any coupling J, the quasi-particle energy 
gap at half-filling is found to be 

A = -7rt + [(7rt)2 + JV4]^/2_ (24) 

The above result indicates that the system is always in in an insulating state for any nonzero 
coupling J and h = 0. Hence, we conclude that \Jc\ = 0+, in spite of the chirality of 
conducting band ( corresponding to long range hopping matrix ). In the strong interaction 



limit, the energy gap is A ?» | J|/2. In the weak interaction hmit, the energy gap becomes 
A fti JV(87rt). 

In summary, the stability of the insulating phase of the ID chiral Kondo lattice model at 
half-filling is studied, within the framework of self-consistent variational theory. It is found 
that the insulating phase is stable for any nonzero impurity-electron interaction, and the 
critical point of metal-insulator phase transition is at J = in zero external magnetic field. 
It would be very interesting to validate the conclusion of the critical point through other 
approaches. Further work is necessary for study of the behavior of the chiral Kondo lattice 
model in nonzero magnetic field. 

We would hke to thank D. Baeriswyl and H. Kunz for conversations. We are also grateful 
to H. Carruzzo for correspondences. This work was supported by the Swiss National Science 
Foundation. 
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